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We present a model-independent description of low-density neutron matter based on the virial 
expansion. The virial equation of state provides a benchmark for all nuclear equations of state 
at densities and temperatures where the interparticle separation is large compared to the thermal 
wavelength. We calculate the second virial coefflcient directly from the nucleon-nucleon scattering 
phase shifts. Our results for the pressure, energy, entropy and the free energy correctly include the 
physics of the large neutron-neutron scattering length. We find that, as in the universal regime, 
thermodynamic properties of neutron matter scale over a wide range of temperatures, but with a 
significantly reduced interaction coefficient compared to the unitary limit. 
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I. INTRODUCTION 

Matter composed of neutrons and protons exhibits 
remarkable properties. Cold neutron matter is close 
to a scaling regime for densities n < 1/10 ng, where 
no = 0.16 fm~^ denotes nuclear saturation density. In 
this regime, thermodynamic properties of neutron mat- 
ter scale only with the density, and for the energy per 
particle E/N one has approximately 0,13 
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where m is the nucleon mass. In contrast, below satura- 
tion density, the energy per particle of symmetric nuclear 
matter is independent of the density E/A « — 16MeV, 
and nucleons cluster into nuclei and larger structures. 
The physics of nuclear matter is therefore a crossover 
from a gas of nucleons to homogeneous matter, where nu- 
clei and larger clusters coexist with the nucleon gas over 
a wide range of intermediate densities. Nuclear cluster- 
ing is due to the competition of nuclear binding, entropy 
and the Coulomb repulsion. 

With this rich interplay, it is very important to de- 
velop a reliable formalism for nucleonic matter over den- 
sities and temperatures that are experimentally or ob- 
servationally accessible. In this paper, we present a 
model-independent description of neutron matter based 
on the virial expansion. We have previously applied the 
virial expansion to low-density nuclear matter composed 
of neutrons, protons and alpha particles The virial 
equation of state presents a controlled application of nu- 
cleon scattering data to densities and temperatures where 
the interparticle separation is large compared to the ther- 
mal wavelength. For neutron matter, the relevant densi- 
ties are n < 4 • 10" (T/MeV)^/^ g/cm^, where T is the 
temperature. The virial equation of state therefore pro- 
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vides important constraints on the physics of the neutri- 
nosphere in supernovae, where one has T « 4MeV and 
n - 10" g/cm^ ~ 1/1000 no 0, Q . 

Properties of dilute Fermi gases can also be studied in 
laboratory experiments with trapped atoms. Cold atom 
experiments provide exciting access to neutron matter 
in the universal low-density regime, where the scattering 
length fls is large compared to the interparticle separa- 
tion fepOs 3> 1 and the effective range is small kpre ^ 1- 
Under these conditions, there are no length scales associ- 
ated with the interaction, and the only dimensionful scale 
is the Fermi momentum fcp. Therefore, the system will 
exhibit universal behavior, where all macroscopic observ- 
ables are given by powers of kp multiplied by universal 
factors. For example, under these conditions the energy 
per particle E/N for cold gases of atomic ^Li, or 
neutrons with equal populations of two spin states is 
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where the universal factor ^ is a number. For neu- 
trons, the scattering length is unnaturally large with 
flriri = —18.5 ± 0.3 fm (for a recent review see f6|). The 
effective range is expected to be approximately charge- 
independent, and thus r„„ = 2.7 fm Consequently, 
the universal regime is restricted to extremely low neu- 
tron densities n = fc|/(37r^) < 10~'*fm"^. 

Dilute Fermi gases with resonant interactions were re- 
alized for the first time by O'Hara et al. in 2002 |^. 
In this and subsequent experiments, the universal fac- 
tor ^ was determined by extracting an equation of state 
from the properties of the atomic cloud |3, ITol [lll |. 
This leads to ^ = 0.51 ± 0.04 ^ ~ 0.7 [13 and ^ = 
0.27 

io 09 E3' temperatures in units of the Fermi tem- 
perature T/Tf « 0.05 (except for T/Tp « 0.6 in 1^). To 
date, the most reliable theoretical results for the univer- 
sal equation of state are from T — 0, fixed- node Green's 
function Monte Carlo simulations, ^ = 0.44 ± 0.01 
and ^ = 0.42 ±0.01 JlJ. 

The virial expansion has been applied to cold atoms in 
the vicinity of Feshbach resonances by Ho et al. 
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The resulting universal equation of state for high tem- 
peratures describes many observed properties of resonant 
Fermi gases, with a universal virial coefficient that is in- 
dependent of temperature. In this paper, we study to 
what extent the virial equation of state of low-density 
neutron matter approaches this unitary limit. 

There are many theoretical approaches to low-density 
neutron matter. Recently, the equation of state of neu- 
tron matter at zero temperature has been calculated us- 
ing an effective field theory for large scattering length 
and large effective range 2]. This provides a model- 
independent description of neutron matter for densities 
where the interparticle separation is comparable to the 
effective range. The effective field theory calculation is 
particularly transparent, since the couplings are fitted 
directly to the scattering length and the effective range. 

In addition, there are microscopic calculations start- 
ing from nucleon-nucleon (NN) and three-nucleon inter- 
actions that reproduce NN scattering and selected few- 
nucleon data. One of these are Hartree-Fock calcula- 
tions Jldl with the model-independent low-momentum 
interaction T^owfc |l7j |. These results are promising, be- 
cause low-momentum interactions offer the possibility of 
a perturbative and thus systematic approach to nuclc- 
onic matter |l8j . Moreover, for neutron matter at subnu- 
clear densities, the Hartree-Fock results agree with com- 
plicated Fermi hyper-netted chain [l^ [20l | and Brueck- 
ner 0,1131 calculations. However, these approaches may 
not be reliable for very low densities, where the large 
scattering length physics is important. This is especially 
problematic for Brueckner calculations, since the hole- 
line expansion relies on Pauli blocking effects being dom- 
inant, so that all large structures are dissolved in the 
medium. Note that the equation of state calculation of 
Buchler and Coon l2^ is close in strategy to the virial 
equation of state, but it takes into account Pauli blocking 
on the phase shifts and thus uses a model NN interaction. 
There are also a number of Skyrme-type 0, |2^ and 
relativistic mean- field [27LI28. 29] parametrizations of the 
nuclear energy functional, which are used to calculate 
ground-state energies and densities of intermediate-mass 
and heavy nuclei. However, the resulting energy func- 
tional for neutron matter is not well constrained from 
fits to finite nuclei. 

Finally, there are promising lattice simulations for neu- 
tron matter using effective field theory, where the cou- 
plings are regularized on the same lattice [s^. These 
calculations are presently limited to small lattices and 
to the lowest-order contact interaction, which is fitted to 
reproduce the scattering length. In this paper, we study 
how the virial equation of state depends on NN scatter- 
ing properties and in particular on the physics beyond 
the large scattering length. As a result, we assess how a 
low-order truncation in the effective field theory impacts 
the equation of state of low-density neutron matter. This 
provides a valuable check for the lattice results. 

The virial equation of state is a general, model- 
independent equation of state for a dilute gas, provided 



the fugacity z = e^^'^ is small. Here denotes the chem- 
ical potential. An additional assumption is that the sys- 
tem is in a gas phase and has undergone no phase transi- 
tion with decreasing temperature or increasing density. 
Under these conditions, the partition function can be 
expanded in powers of the fugacity. The second virial 
coefficient 62 describes the term in this expansion 
and is directly related to the two-body scattering phase 
shifts 1^ I33. Little is known about the third virial 
coefficient, which describes the term |33,I3- We em- 
phasize that the virial expansion is not a perturbative 
expansion in powers of k-pds- A great advantage of the 
virial formalism is that it includes both bound states and 
scattering resonances on an equal footing. It correctly 
predicts that thermodynamic quantities, such as the pres- 
sure, are continuous as the interaction is changed to con- 
vert a low-energy scattering resonance into a weakly- 
bound state. This continuity has been demonstrated ex- 
perimentally with cold atoms in the crossover region of 
a Feshbach resonance (see e.g., ^3) ^"^^ ^-Iso shown 
theoretically by Ho and Mueller using the viral expansion 
across Os = ±00 [T3 |. 

This paper is organized as follows. We briefly intro- 
duce the virial equation of state in Section ^] Further 
details can be found in 0. In Section IhII we present 
results for the second virial coefficient, the pressure, en- 
ergy, entropy and the free energy of low-density neutron 
matter. We also study the dependence of the virial coeffi- 
cient on NN scattering properties and show that neutron 
matter scales to a very good approximation. Finally, we 
conclude in Section Hvl 

II. VIRIAL EQUATION OF STATE 

For the virial equation of state we expand the grand- 
canonical partition function or the pressure in a power 
series of the fugacity 

2T 

P = — (z + z25„ + ^3^(3) ^ 0(^4)) ^ 

where A denotes the nucleon thermal wavelength A = 
{2iT /mTY^'^ and 6„, bn^ are the second and third virial 

coefficients for neutron matter respectively. We will in- 

(3) 

elude bii only to make an error estimate for our virial 
results. The density follows from differentiating the pres- 
sure with respect to the fugacity n — z/T {dzP)v,T find 
is given by 

?i = ^ (z + 2z26„ + Sz^fcCf ) + 0(^4)) 

Therefore, the fugacity expansion is an expansion in pow- 
ers of nA"^, and for finite temperatures, it is valid to much 
higher densities than the k-pag expansion. In this work, 
we truncate the virial expansion after second order in 
the fugacity. This leads to an equation of state that is 
thermodynamically consistent. 
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The dependence of the density on z can be inverted. 
This gives the virial equation of state directly in terms of 
density and temperature P = P(^z{n, T), T). In practice, 
we directly, without this inversion, generate the equa- 
tion of state in tabular form for a range of fugacity val- 
ues. This maintains the thermodynamic consistency of 
the virial equation of state. 

The second virial coefficient is related to the partition 
function of the two-particle system X^statcs e~^^/"^, where 
the sum is over all two-particle states of energy E2 ■ This 
sum can be converted to an integral over relative momen- 
tum k weighted by the density of states of the interacting 
two-particle system. The density of states, and therefore 
the second virial coefficient, can then be expressed in 
terms of the scattering phase shift 6{k), summed over all 
allowed partial waves iSlj, |32| . 

For neutron matter, the second virial coefficient 6„ is 
given by 



bn{T) 
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where — 2~^/^ is the free Fermi gas contribution and 
is the sum of the isospin- triplet elastic scatter- 
ing phase shifts at laboratory energy E. This sum is 
over all partial waves with two-particle spin 5, angular 
momentum L and total angular momentum J allowed by 
spin statistics, and includes a degeneracy factor (2J-I- 1), 
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Note that we have neglected the effects of the mixing 
parameters due to the tensor force. We expect that their 
contributions to the second virial coefficient describing 
spin-averaged observables vanish. 

Finally, the entropy S and the energy E are obtained 
from the virial equation of state using thermodynam- 
ics 0,122. The entropy density s = S/V follows from 
differentiating the pressure with respect to the tempera- 
ture s = {dTP)fj.- This leads to 



— -nlogz + —z 



(7) 



where b'^iiT) = dbn{T)/dT denotes the temperature 
derivative of the virial coefficient. The energy density 
e = E/V can be calculated from the entropy density. 



e — Ts + n/i — P 



lp-^^z\' 
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For completeness, the entropy per particle S/N , energy 
per particle E /N and the free energy per particle F/N 
are given by 



S _ s E _ e F _ f 

with the free energy density f — e ~ Ts. 



(8) 
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FIG. 1: (Color online) The total phase shift S*°\E) ver- 
sus laboratory energy E obtained from including successively 
higher partial waves. 



III. RESULTS 
A. Virial Coefficients 

We first calculate the second virial coefficient 6„. We 
use the T = \ np phase shifts obtained from the Nijmegen 
partial wave analysis j35j] . This neglects the small charge 
dependences in nuclear interactions. We have included 
all partial waves with L ^ 6 and the resulting total total 
phase shift (5*°' {E) is shown in Fig. The total phase 
shift displays both the prominent low-energy '^Sq reso- 
nance and the significance of higher partial waves. We 
observe that (5*°' [E) is very weakly dependent on energy 
with (5*°' « 65 degrees over the entire range E ^ 350 MeV 
of the Nijmegen partial wave analysis. As shown in Fig.^ 
the decrease of the ^Sq phase shift due to the large ef- 
fective range is compensated by the contributions from 
higher angular momenta. For reference, in the unitary 
limit of a S-wave resonance (a^ = ioo), the phase shift 
is energy- independent with 61^^ = 90 degrees. 

In the unitary limit, the second virial coefficient is inde- 
pendent of the temperature and given by 62 = 3/2^/^ = 
0.53 For neutron matter, our results for the virial co- 
efficient hn and its derivative Tb'^{T) are listed in TableQ] 
In addition to the full results, we also give the virial co- 
efficient obtained only from the S-wave scattering length 
a„p and including the effective range contribution. As 
expected, 6„ is dominated by the large S-wave scattering 
length, but effective range and higher partial wave contri- 
butions are significant even for these low temperatures. 
As a result of the weak energy dependence of 5^°^{E), we 
find that the second virial coefficient is approximately 
independent of temperature over a wide range, and con- 
sequently Th'^ « 0. The value we obtain for 6„ w 0.31 is 
40% reduced compared to the universal value 62 = 0.53. 
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TABLE I: The second virial coefficient b„ for different temperatures. In addition to the fuU resuhs, we also give results for the 
virial coefBcient calculated in the unitary limit and from including only the S-wave scattering length a„p, with and without the 
effective range re contribution (a„p = — 23.768fm and re = 2.68 fm 7]). The results labeled CIB take into account the effects 
due to charge-independence breaking (CIB) on the scattering length with a„„ = — 18.5fm. 
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TABLE II: The second virial coefficient 6„ for temperatures 
T > 25MeV. As discussed in the text, we estimate an error 
of < 5% for these higher temperatures due to the truncation 
of the integration over the phase shifts at _B ^ 350 MeV. 
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In Table we also study the effects of charge- 
independence breaking (CIB) on the scattering length. 
Due to the lack of neutron-neutron scattering data, we 
estimate CIB effects by subtracting the virial coefficient 
calculated only with a„p from the full 6„, and then add 
the virial coefficient obtained from the neutron-neutron 
scattering length a„„ = — 18.5fm. We find that CIB 
effects are largest for T < 5 MeV and lead to a 10% re- 
duction of the second virial coefficient. 

Finally, we extend our results for 6„ to temperatures 
T ^ 25 MeV in Tabled We give these results separately, 
because there is a small error for these higher tempera- 
tures due to the truncation of the integration over the 
phase shifts at i? < 350 MeV (the extent of the partial 
wave analysis). Assuming the total phase shift is con- 
stant, we vary the energy cutoff to > 350 MeV and 
estimate the error to be < 5% up to the highest temper- 
atures given in Table UTI 
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FIG. 2; (Color online) The pressure P versus density n for 
T = 10, 16 and 20 MeV. The dotted error bands for the virial 
equation of state are based on an estimate of a neglected 
third virial coefficient b^^ — ±6„/2. Note that a negative 
increases the pressure. Also shown are the FHNC results 
of Friedman and Pandharipande (crosses) lig|l . The circles 
indicate where the fugacity is z = 0.5, and the error bands 
are shown for 2: ^ 1. For this density range, the fugacities are 
z < 1.65, 1.01 and 0.79 for T = 10, 16 and 20 MeV. 



B. Pressure 

Our virial results for the pressure are shown in Fig. [2| 
for temperatures T = 10, 16 and 20 MeV. The third 
virial coefficient can be used to make a simple error es- 
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FIG. 3: (Color online) The pressure P versus density n calcu- 
lated for r = 1, 4, 8 and 20 MeV and scaled to a temperature 
of T' = 4 MeV (for details see text) . The lower curves are re- 
sults for the scaled pressures obtained with a virial coefficient 
calculated only from the S-wave scattering length a„p. The 
latter curves have been shifted for better readability. The cir- 
cles indicate where the fugacity is 2 = 0.5, and for this density 
range, the fugacity is z < 0.87. 



FIG. 4: (Color online) The energy per particle E/N versus 
density n for T = 4, 8 and 20 MeV. In order to clearly show 
the interaction effects we have subtracted the free kinetic en- 
ergy 3T/2 from all curves. The energy per particle obtained 
using the universal relation e — 3P/2 with the corresponding 
virial pressures are indistinguishable from these curves. The 
circles indicate where the fugacity is z — 0.5, and for this 
density range, the fugacities are z < 1.60, 0.77 and 0.25 for 
T = 4, 8 and 20 MeV. 



timate of neglected terms in the virial expansion. The 
Pauli principle prevents three neutrons to interact in the 
S-wave and there is no three-neutron bound state. More- 
over, 3N forces are very small in low-density neutron 
matter. Therefore, we consider a third virial coefficient 
\bn \ < bn/2 reasonable. Note that for the ideal Fermi 
gas bi^-* — 3^^/^ — 0.06. In Fig. [3 we show the resulting 
error bands, which are small even up to z ^ 1. We also 
compare our results to the the microscopic Fermi hyper- 
netted-chain (FHNC) equation of state of Friedman and 
Pandharipande We find that the FHNC results 

are in very good agreement with the model-independent 
virial equation of state. Moreover, it is intriguing that 
the FHNC results lie within the error band for the virial 
pressure even for fugacities z > 0.5. In contrast, we have 
found a disagreement of the FHNC results with the virial 
equation of state for low-density nuclear matter due to 
clustering Q. 

C. Scaling and Energy 

If the virial coefficients are temperature independent, 
then the power series in the fugacity for the pressure and 
the density will have no explicit temperature dependence. 
Consequently, the dependence of the pressure on density 
and temperature will be given by 

P{n,T) = T'^/^ f{n/T^/^), (9) 

where f{x) is a scaling function oin/T^/'^. If this scaling 
relation holds, one can predict the pressure T') at 



a new temperature T' from P(n, T) through 

P{n',T')={^y\{n,T), (10) 

with n' = (T7T)3/2 n. In Fig. El we demonstrate that 
low-density neutron matter scales according to Eq. H10|l . 
The predicted pressures for T' = 4 MeV obtained from 
virial pressures calculated for T = 1,8 and 20 MeV are in 
excellent agreement with the unsealed T — A MeV virial 
pressure. Fig. |21 also gives results for the scaled pressures 
(T 1, 8 and 20 MeV scaled to T' = 4 MeV and unsealed 
T — A MeV) obtained from the virial equation of state 
with bn calculated only from the S-wave scattering length 
a„p. Although Tabled shows that 5„ increases by 35% 
from r = 1 to 20 MeV in this case, the scaling symmetry 
continues to hold to a good approximation. 

In the scaling regime, 6^ = and therefore the energy 
density, Eq. (jS)), is given by 

e = 3P/2. (11) 

This is a general thermodynamic relation in the universal 
regime |36| . In Fig. 01 we show the energy per particle 
for T = 4, 8 and 20 MeV. In order to separate the inter- 
action effects we have subtracted the free kinetic energy 
3T/2 from all curves. We have also calculated the energy 
per particle from the virial pressures using the universal 
relation Eq. (|ll|l . The resulting curves for the energy per 
particle are indistinguishable from those shown in Fig. 0] 
This demonstrates that the b'^ term is indeed negligible. 
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FIG. 5: (Color online) The entropy per particle S/N versus 
density n for T — 10, 16 and 20 MeV. Also shown are the 
FHNC results of Friedman and Pandharipande (crosses) [l9| . 
The circles indicate where the fugacity is z = 0.5. 



FIG. 6: (Color online) The free energy per particle F/N ver- 
sus density n for T — 10, 16 and 20 MeV. Also shown are the 
FHNC results of Friedman and Pandharipande (crosses) [l9|. 
The circles indicate where the fugacity is z = 0.5. 



While neutron matter scales as in the unitary limit, the 
thermodynamic properties depend on the value of the 
virial coefficient, and the latter depends on the physics 
of NN scattering. 

Finally, we emphasize that the universal relation for 
the energy density, Eq. (fTT)) . holds independent of the 
value of the universal interaction coefficient given by ^ 
for T — and 6„ for higher temperatures. In fact, 
Eq. is valid for T = 0, where the energy per par- 
ticle only scales with density, E/N ~ ^ n^^^ according 
to Eq. 101, and for the virial regime, where the en- 
ergy per particle scales with density and temperature 
E/N ~ ^(T/n2/3) ^2/3 ^ l{T/n'^/^)T, with Fermi tem- 
perature Tp ~ n^/^. The latter scaling is very explicit to 
lowest order in the density, where the energy per particle 
is given by 

|«^r(l-6„AV2). (12) 



D. Entropy and Free Energy 

We present our virial results for the entropy per par- 
ticle in Fig. Elfor temperatures T = 10, 16 and 20 MeV. 
The FHNC results are again in very good agreement with 
the model-independent virial equation of state, with the 
exception of two FHNC points that seem to either reflect 
numerical instabilities or are typos in jl9| . As indicated 
by the circles for z = 0.5 in Fig. [3 we observe that the 
range of validity of the virial equation of state is bounded 
by a constant entropy per particle, independent of the 
temperature. This is easily understood from the expres- 
sion for the entropy per particle in the scaling regime 



with 6; = 0, 

N 2\ l + 2zbnj ^ ^ ' 

Therefore, if one takes as the range of validity of the virial 
equation of state z < 0.5, this is equivalent to S/N > 2.9 
for bn = 0.31, in agreement with Fig. [S] For complete- 
ness, we also show the free energy per particle in Fig. 
As for the pressure and the entropy, the FHNC results 
agree well with our model-independent virial predictions. 

IV. CONCLUSIONS 

We have presented the virial equation of state of low- 
density neutron matter. The virial expansion provides 
a systematic way to include strong interactions with a 
well-defined range of validity. The resulting virial equa- 
tion of state has a simple parametric form and is thermo- 
dynamically consistent. We have calculated the second 
virial coefficient directly from the NN scattering phase 
shifts. Therefore, the virial equation of state is model- 
independent and sets a benchmark for all nuclear equa- 
tions of state at low densities. 

We have found that the second virial coefficient 6„ is 
approximately temperature- independent. This is the re- 
sult of an intriguing cancellation. The decrease of the 
^So phase shift due to the effective range is compensated 
by the increase of the higher angular momentum phase 
shifts. Therefore, while 5„ is generally dominated by the 
low-energy ^Sq resonance, contributions from higher par- 
tial waves are significant even for these low temperatures. 
The virial coefficient 6„ « 0.31 is reduced by 40% com- 
pared to the unitary limit, where 62 = 0.53. 
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The virial equation of state was used to make model- 
independent predictions for the pressure, energy, entropy 
and the free energy of low-density neutron matter over 
a wide range of densities and temperatures. Our results 
include the physics of the large neutron scattering length 
in a tractable way. The range of validity of the virial 
equation of state is given by the fugacity. With z < 0.5, 
one has for the entropy per particle S/N > 2.9 or for the 
density n < 4-10" (T/MeV)3/2 g/ cm . The virial expan- 
sion thus provides important constraints on the physics 
of the neutrinosphere in supernovae. 

We have made simple error estimates of the virial equa- 
tion of state by studying the effects of a neglected third 
virial coefficient \bn \ ^ bn/2. The resulting error bands 
are small. For a better error estimate, it is important to 
have a reliable calculation of the third virial coefhcient. 
The FHNC results of Friedman and Pandharipande 
are in very good agreement with the model-independent 
virial equation of state. This is in contrast to our find- 
ings for low-density nuclear matter Q. The FHNC and 



virial results agree within the estimated errors even for 
densities, where z > 0.5. 

The temperature independence of &„ leads to a scal- 
ing symmetry of low-density neutron matter. As in the 
universal regime, thermodynamic properties of neutron 
matter scale over a wide range of temperatures. This 
extends the approximate scaling of the T = equation 
of state according to Eq. to finite temperatures. Fi- 
nally, model-independent predictions for the density and 
spin response of low-density neutron matter and a de- 
tailed comparison with nuclear lattice calculations will 
be discussed in future work. 
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